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<N : 1 ABSTRACT 

^ ■ A gauge theory of pions interacting with rho-mesons at elevated temperatures is used to calculate 
the pressure in a hot pion gas. No reference is made to the pion's status as a QCD Goldstone 

I boson. The role of the pion is merely that of a carrier of an SU(2) symmetry, gauged to create a 

^Sl ' vector-meson interaction, the rho playing the role of the interacting vector particle. The results are 

. in rough agreement with much more elaborate calculations, both of the purely hadronic variety, 

I and those that invoke quark-gluon degrees of freedom. The quark-gluon and purely hadronic 

' calculations seemingly lead to very similar predictions which are in accord with receent data from 

■ RHIC. The results motivate the question as to whether the two descriptions are dual to each other 

\^ I in the sense of being alternate models, each sufficient to explain the observed data. 

o ■ 

Experimental groups at RHIC have recently summarized their results from three years of measure- 
r~| ! ments on gold-gold collisions with center-of-mass energies per nucleon pair of 130 and 200 GeV. 

I I shall focus, for definiteness, on results from the 200 GeV per nucleon pair runs reported by the 

PHENIX experiment p. 

Heavy ion collisions, even before the RHIC experiments, were found to produce clouds of hun- 
I dreds of particles, with a predominance of pions, per collision l^. One can imagine that, in 

accordance with Bjorken's description [S], the colliding ions pass through each other like porous 
pancakes and vanish down the beampipe, leaving behind a fireball of pions and other hadrons. 
Like-pion correlation measurements at RHIC reported in reference pOj (Section 3.4) suggest that 
the fireball has about the size of a gold nucleus (radius ~ 6 fm. There are about 1000 pions per 
unit rapidity produced in a central collision of gold nuclei with ^/sNN = 200 GeV at RHIC jl], 
leading to an estimated pion density in the fireball of about 2 fm~'^. These fireballs are a new 
form of matter (characterized by Shuryak |5j as "the pion liquid" - partially in recognition of the 
high density of particles in the clouds). 

A task for the theorist, faced with the ephemeral presence of this rare form of matter, is to 
obtain its equation of state and investigate its observable properties. A number of authors have 
attacked various aspects of this task during the past several decades, see e.g. [3], 0, |HI, jHl; 
|10j . and papers cited therein. The common approach is through the relativistic thermodynamics of 
an interacting pion gas. Reference 9 is based upon a hadronic effective Lagrangian for the pionic 
interactions, with in-medium modification of the single pion propagator. Reference |10j bases 
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its vrvr interaction on chiral perturbation theory. Both approaches reproduce experimental (zero- 
temperature) scattering data so that thermodynamic variables calculated in the two approaches 
are expected to, and appear to, agree by virtue of a theorem of Dashen, Ma and Bernstein 

Reference ^0] includes density considerations by noting that pi-pi scattering is almost purely 
elastic in the energy range of interest, so that pion number can be treated as a conserved quantity. 
This observation is used to justify the inclusion of a chemical potential corresponding to pion 
number. See also, Welke, et al. ^ 

The approach in the present paper exploits a remark made long ago in reference jl2j . The 
authors say of their chiral perturbation theory results for pi-pi scattering (p. 192): "We show that 
in a systematic low-energy expansion in powers of the momenta the p does not play any special role 
- its presence only manifests itself in the values of the low energy constants" . But we learn from 
reference that the knowledge of the S-matrix is sufficient for calculating the grand canonical 
potential of any gaseous system. I take this as meaning that a thermal field theory of pi's and rho's 
should provide a "pretty good" description of a hot pion cloud. The meaning of "pretty good" may 
be gleaned from the fact that p-exchange in the tttt- scattering t-channel gives respective I-spin 
and 2 scattering lengths of .36m~^ and — .18m~^, whereas a recent precision determination of the 
lifetime of the pionium atom by the Dirac Collaboration ^S] gave |ao — 02! = .264^*^ 020 

3 The Lagrangian 

As proposed in 1967 by Lee and Nieh and Weinberg ^3], I start by providing the pion field 
with a local SU(2) gauge transformation to generate a gauge-invariant coupling to an isovector 
vector field. I then break the SU(2) symmetry by providing the vector field with the mass M of 
the p - meson. The resultant broken symmetry is sometimes referred to as a "hidden local gauge 
symmetry" |lf)| . The "hidden symmetry" manifests itself by preserving coupling-constant relations, 
despite the breaking. 

Introducing the vector field and the covariant derivative 



(1) 



where 



T xT = iT 



(2) 



the pion Lagrangian density is 

jCtt = ~ 2^*^^ ' ~ m^vf • 7?) 
where the vector symbols refer to SU(2) space. The total Lagrangian density is then 



(3) 



J^TOTAL = - \Ff,, ■ F^'^ - hl\V^ ■ y^) 



(4) 



with 



V 



(5) 



2 



Figure 1: Feynman diagrams for p-polarization tensor; p- full line, pion- dotted line 



The interaction Lagrangian density is the part of Ctotal that is proportional to g and g^. 

CiNT = -gV^ ■ [v? X d^vf] - ^<72y^ • [vr X {V^ X vf)] (6) 

4 THERMAL MASS AND WIDTH OF THE RHO 

The order g'^ contribution from the interaction described by Eq. El to the thermal p-polarization 
tensor can be calculated using the methods described by Kapusta ^7], Le Bellac ^H]) and Pisarski 
jl9j The relevant Feynman diagrams are shown in Fig. 1. I evaluate them in the Matsubara 
formalism |2Uj . 



Fig.l Feynman diagrams for the second order p-polarization tensor. 
The p-propagator in Euclidean momentum space is: 

Diq)t.u = L{q)f,^-^^j^ (7) 

with 

m,. = V - ^ (8) 

The choice of Landau gauge for the propagator insures that the p will have three components both 
on- and off-shell. 

The squared mass of the p at temperature T, including one-loop order pion contributions, is 
given by: 

M{pf-M'^ = -^^T^y{fd^kd^k'A{k)A{k')f^dy[e'y^P'''''''~'''o^d{p-k 
3{2iTy J Jo 

n,n 

p-^ 2, 6p'^ 

(9) 

where the argument "0" in a propagator refers to the spacetime origin, A(A;) is (inconsistently) 
the momentum-space pion propagator, and D^^'^{p) is the momentum-space p-propagator. The 
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Figure 2: Mass and width of the p vs temperature 



indices n, n' refer respectively to the discrete Matsubara [201 modes of the fourth components of 
the Euchdean four- vectors "A:. A;'" 

Further details of the p-mass calculation are provided in Appendix A. The temperature depen- 
dence of the mass and width of p -mesons at rest in a heat bath are repectively shown in Fig. 
2. 



It is evident from Eq. ^Jthat the effect of the heat bath is a negative contribution to the mass of 
the p which turns out to be small at temperatures below 200 MeV. The width is Bose-enhanced, 
as previously discussed in [2J, in accordance with Eq. 1171 These results are in at least qualitative 
agreement with results from the NA60 experiment at the CERN SPS [22] 
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5 EQUATION OF STATE OF A HOT PION GAS 



The pressure p of a hot pion gas, to second order in the ir — p couphng, may be written in the form 

(im, d) 

l^ = A, + 39^A (10) 
The free-pion contribution Aq is given by 

^o = ^fniPu;)^ (11) 

with uj denoting the energy of a pion of momentum k, (3 the inverse of the temperature T, and 
n{x) the Bose-Einstein distribution function (e^ — 1)~^. 

The order pion-contribution to A in Eq. can be divided into two contributions, as shown 
in Fig. 1. The corresponding expression may be written 



(fkd^Pj 



A = S / ^—f[Aip,k) + B{p,k)] (12) 



with 



A{p,k) = ^ J d^k'6{k + k' +p) j\Te'^Po)+ko+K)^^ 
A{k)A{k')Df'''{p){k - k')^{k - k'U 



(13) 



and 



S(p,fc)=r2Z)-(0)A(0) (14) 

The summation is over the respective Matsubara modes of the fourth components oi p,k,k' , and I 
have symmetrized in the momenta k and k' . Some details of the calculation are in Appendix B. A 
graph of the resulting pressure-temperature relation is shown in Fig. 3, where it is compared with 
the corresponding results obtained by Rapp and Wambach. [^. Fig. 3 shows that the simple model 
presented here overestimates the pressure by about 30% at a temperature of 170 MeV compared 
with the model in reference [S] 



The result obtained in ref. ^ is almost identical with a QCD-related calculation by Nikonov, 
et al. |23j . That calculation treats the "fireball" aftermath of the heavy-ion collision as a "mixed- 
phase" combination consisting of both a quark-gluon phase and a hadronic phase. 

6 The Pion Gas at RHIC 

The aftermath of a heavy ion collision is an expanding fireball of pions and other debris. Transverse- 
momentum spectra of it's, K's p's and p's are plotted and tabulated in ref. ^ for ^/sNN = 200 
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T(MeV) 



Figure 3: Pressure v Temperature in a pion gas; dashed line-this calculation, dotted line- Rapp & 
Wambach 

GeV gold on gold collisions at RHIC. The transverse momentum of a debris particle is compounded 
from the expansion of the cloud and the thermal motion of the particles in the cloud. Siemens and 
Rasmussen |23] pointed out three decades ago that the heavy particles in the expanding cloud 
can be used to measure the expansion velocity of the cloud. Their point was that the peak of a 
heavy-particle spectrum represents the momenta of those particles with zero transverse thermal 
velocity. 

Fig. 4 is a plot of the transverse-momentum spectrum of mesons taken from Table XVII 
of ref. 15 (centrality 0-5%, error bars omitted), 



Although the peak is not included in the reported measurement, we can roughly estimate that it 
would be near, or a bit below, .35 GeV. Taking the charged-K mass as .5 GeV, results in a /32 of 
.32, or a bit less, for the squared speed (divided by c) of the expanding cloud. This result is close 
to the limit of 1/3 for a gas of highly relativistic particles and about equal to the expansion 
velocity of a free-pion gas at a temperature in excess of about 150 MeV, as reported by Welke, et 
al. 0. 

The observed slope of the kaon spectrum corresponds to an apparent temperature of Tapp ~ 290 
MeV. Correcting this value by transforming to the frame of the expanding gas cloud result |2l] 
results in a Kaon temperature T = \Jj^Tapp ~ 150 MeV, near the range of temperatures estimated 
in ref. pp. 
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Figure 4: Transverse momentum spectrum of mesons from Au+Au collisions at ^/sJfN GeV at 
RHIC 

7 CONCLUSIONS 

The implicit questions underlying the present work are: What are the limits of the theorem of 
Dashen, et al. (llj . that makes the equation of state of a collecttion of hadrons an offspring of 
a hadronic S-matrix? Can purely hadronic descriptions of the thermodynamics of the "fireball" 
created by heavy ion collisions account, at least qualitatively, for the results observed at RHIC? 
PP? Or, must quark-gluon degrees of freedom be invoked to properly describe the RHIC results? 

The questions are posed in the context of a simple model, a gauge theory of pions interacting 
with /3-mesons at elevated temperatures. No reference is made to the pion's status as a QCD 
Goldstone boson, in contradistinction to the works cited in the Introduction. The role of the pion 
is merely that of a carrier of an SU(2) symmetry, gauged to create a vector-meson interaction, the 
p playing the role of the interacting vector particle. The results are in rough agreement with much 
more elaborate calculations, both of the purely hadronic variety, and those that invoke quark-gluon 
degrees of freedom. The quark-gluon and purely hadronic calculations seemingly leading to very 
similar predictions motivates the question as to whether the two descriptions are dual to each other 
in the sense of being alternate models, each sufficient to explain the observed data. There are other 
indications in the published literature that would seem to motivate the same question. 

Kapusta, Lichard, and Siebert |2ni calculate the photon emission from a hot hadronic matter 
and from a quark-gluon plasma at temperatures from ^ GeV to several GeV and find them to be 
the same. Scherer and Fearing [27] show that different off-shell form factors for s- and u- channel 
pole diagram contributions to pion compton scattering, calculated in chiral perturbation theory, 
lead to the same S-matrix, in accordance with Haag's Theorem j28j . 

Haag proved that in a theory with an S-Matrix (which excludes both QED and QCD) there 
is no unique off-mass-shell extension. The necessary conditions for this lack of uniqueness are 
unknown, to the best of my knowledge. It may be that modern high-energy experiments, which 
count jets, rather than isolated free particles, avoid Haag's theorem and admit the presence of 
unique Lagrangians. Whether the study of heavy ion collisions can help explore the reach of Haag's 
theorem remains to be seen. 
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A RHO-MASS DETAILS 



Summation over the Matsubara modes and evaluation of the angular integrals in Eql^then gives, 
for p-mesons at rest 

,2 



M{p = Of = m2--^/ edk 

12-11^ Jo 



Jo uj'^ sinh{(3uj / ly 

(6 - ^)coth{(5E/2) 3coth{l3uj/2) kcoth{kp/2)) 
2E 4^ 2M2 ^ 

(15) 

E and to in EqElare the respective energies of a p or pion having spatial momentum of magnitude 
k. The y integration in the preceding equation results in a term ^'^^°*^(^^^/^) ^ xhe coth(x) functions 

in the equation contain both zero temperature parts and temperature dependent parts according to 
the identity cot/i(|) = 1 + 2n(x), where n(x) is the Bose-Einstein distribution function (e^ — 1)^^. 
The real zero-temperature part contributes to the infinite mass renormalization and is incorporated 
into the constant term M^, (see e.g. |17j . p. 40). With the understanding that the constant term 
refers to the renormalized mass, then (and noting that the last term in Eq. EJcan be evaluated in 
closed form), the expression to be evaluated numerically becomes 

M(p = of = m2 - At / k'^dk 







2kM^(5) ^ i6-^)n{(3E) ^ 3n(/3u;) ^ ^ vr^T^ 



CO UJ 



2+p2/4] 2E 4uj ' 15M2 



(16) 



I obtain the width by continuing the fourth component of the p momentum four-vector to a 
positive imaginary plus a real infinitesimal, contributing vr times a delta function multiplying the 
first term of the integrand in Ea llSl The result is 

MT = ImM" = ^^M\l - 4(^)')^[1 + 2n(M/3)]. (17) 

The zero-temperature width given by Eq. El^-grees with that given by Nishijima j29j and determines 
the value of the gauge coupling constant g. The temperature dependent part of the equation predicts 
that the p width broadens with temperature when the p is immersed in a heat bath due to the 
increased density of final states. 



B EQUATION OF STATE CALCULATION 

The p-propagator D{pY^ in Eq El has two poles on the real po-axis, at \p\ and at E{p). It is 
convenient to separate the poles, so that the propagator D{pY'^ becomes -^[^ — pi^j^fA P^^ip)^^ ■ 
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After doing the Matsubara sums and the the angular integrals the resulting equation for the second 

p2 
J'4 ) 



order (in g) contribution to the dimensionless pressure, S^, becomes 



tt^^Att' Jo uj{k) ' ' E{k) 

where 



(19) 



fracn[pioiq)]uj{q) — ^ In i? - kq] + 



n[/jii;(,)] [(M^ - V) q^ 
E{q) ^ 8 ^'^^+'^9^1 + M2J 

(20) 



and 



-?,kn[l3E{k)] (^^^_^^n[/3a;(g)] , ,n[/3i?(g)] . 



The two arguments of the log functions are R = -t^^]l%-%X%+% ^ = (S'-V)E(gj't(g^2fcj' ■ 
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